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We calculate the nucleon contribution to the photon self-energy in a plasma, including the effect 
of the anomalous magnetic moment of the nucleons. General formulas for the transverse and longi- 
tudinal components of the self-energy are obtained and we give explicit results in various limits of 
physical interest. The formulas are relevant for the study of the photon dispersion relations and the 
dynamical susceptibility in a nuclear medium such as the core of a supernova, and has implications 
with regard to the recent suggestion that the Cerenkov process v — > 1/7 can take place in such a 
system. 



I. INTRODUCTION 

The methods of Finite Temperature Field Theory (FTFT) Q have proven to be very useful in the study of the 
matter effects on the propagation of the elementary particles. In these studies, a central role is played by the inverse 
propagator, or cquivalently, the self-energy of the particle. The latter quantity enters directly in the linear part of the 
equations of motion for the effective classical field, and therefore determines the modes that can propagate through 
the medium. In a modern context, this point of view has been stimulated by the influential work of Weldon ^-|}, 
who emphasized the convenience of covariant, real-time, calculations. His approach has been applied to investigate 
the effects of a medium on the properties of photons, gauge fields, and neutrinos, in a variety of physical conditions. 

The case of photons occupies an important place within these type of problems. Its applications range from the 
calculation of cross sections for particle emission by stars, to the determination of the optical properties of solids. 
Under certain conditions of the matter background, the photon dispersion relation in a plasma can be obtained by a 
semiclassical approach based on the Maxwell-Boltzmann-Vlasov equation || . It seems that the first complete quatum- 
field-theoretic treatment of the same problem is due to Tsytovich || . He generalized earlier work along similar lines 
0, obtaining results that are valid under more general conditions than those based on the semiclassical approach. 
The work of Weldon J2| demonstrated that the real-time formulation of FTFT is well suited to study this physical 
system in an efficient and transparent way, and can be extended to, for example, non-abelian gauge fields. 

For photons that propagate in a hot plasma, the medium typically consists of an electron gas superimposed over a 
static positive background. It can be shown that the results of the semiclassical calculation are recovered from those 
of the ficld-thcoretic treatment, by taking the limit in which the photon energy and momentum are small relative to 
the electron mass. Ordinarily, only the electrons play an important role, and the dynamical effect of the nucleons in 
the background are ignored. This is justified within the context of a semiclassical approach. In fact, the contributions 
of the ions and electrons to the dielectric constant have similar form. Since the plasma frequency associated with a 
certain particle is inversely proportional to its mass, the contribution of the ions turns out to be negligible, at least 
in the most common situation where all species share a common temperature |8| . 

In a quantum-theoretic approach, however, the electromagnetic couplings of the nucleons must include their anoma- 
lous magnetic moments. For the case in which the nucleons constitute a degenerate non-relativistic gas, a standard 
treatment yields the Pauli formula for the paramagnetic spin susceptibility, in the limit of zero photon frequency 
(static limit) j|. Nevertheless, a general treatment of the effect in the photon dispersion relation due to the anoma- 
lous electromagmetic couplings of the nucleons seems to be lacking in the literature. Such a study has applications in 
some astrophysical contexts, where systems with dense nuclear matter need to be considered. For example, it has been 
recently pointed out by Mohanty and Samal JlOj ] that, in a supernova core, the photon dispersion relation receives a 
contribution from the nucleon magnetic moment which is of the opposite sign to that of the usual electron plasma 
effect. Under such circumstances, the emission of Cerenkov radiation by neutrinos becomes possible, which can have 
important implications for the energetics of the system. More recently, Raffclt [[ll] has noted that the conclusions of 
Ref. |To| ] are based on the adoption of the static paramagnetic susceptibility formula mentioned above, in order to 
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determine the effect of the nucleons on the photon dispersion relation. However, for a proper account of this effect 
what is required is the full nucleon contribution to photon self-energy (or, equivalently, to the dielectric constant or 
the paramagnetic susceptibility) for any value of the photon momentum and not just in the zero-frequency limit. The 
necessary calculations along these lines are not in the literature, to the best of our knowledge. 

Motivated by the above considerations, in the present work we calculate the photon self-energy in a matter back- 
ground including the contribution from the nucleons, and in particular the effect of their anomalous electromagnetic 
couplings. The calculation is based on the application of FTFT to this problem and complements the existing cal- 
culations of the same quantity for the electron gas alone. The results for the photon self-energy can be equivalently 
interpreted in terms of the dielectric constant and the susceptibility of the system, and in that way we show that the 
well known textbooks results for these quantities are reproduced when the appropriate limits are taken. On the other 
hand, the results we obtain are valid for general conditions of the nucleon gas, whether it is degenerate or not, and 
also for general values of the photon momentum and not just for some particular limiting values. Therefore, they 
are relevant and useful in the study of electromagnetic processes ocurring within dense nuclear matter, such as the 
neutrino Cerenkov radiation in a supernova core mentioned above. 

In Section |n] we give the general 1-loop formulas for the generic contribution of a fcrmion to the photon self-energy, 
including the effect of a possible magnetic moment coupling of the fermion to the electromagnetic field. There it is 
shown that the contribution from any given fermion can be written in terms of just three independent functions, which 
are expressed as integrals over the momentum distribution functions of the fermion. In Section III explicit formulas 
are given for various limiting cases of physical interest. The formulas so obtained are used in Section [V to study the 
photon dispersion relations in a background composed of electrons and nucleons under various conditions, including 
those corresponding to a degenerate or a classical nucleon background. Although the focus of our calculations is 
the real part of the self-energy, there we also summarize the results of the calculation of the imaginary part, which 
determines the damping of the propagating modes. Some details of that calculation are given in Appendix [b| Section 
^ contains our conclusions, and in Appendix ^ the relation between the self-energy and the dielectric function of 
macroscopic electrodynamics is briefly reviewed. 



II. PHOTON SELF-ENERGY 

As shown in Ref. [0, the equation of motion for the effective field in the medium, from which the dispersion 
relations of the propagating modes are determined, is given by 



A v = 



Jim- 



(2.1) 



where 



9\iv — 9fj,u 



(2.2) 



and 



r (eff) 



■ u) + Tr*6(-q ■ u) . 



(2.3) 



In Eq. (2.3) 6 is the step function, stands for the velocity four- vector of the medium and n^u is determined in 
terms of the elements of the 2x2 self-energy matrix in the form 



Re 717^ = Re7rii M „ , 



(2.4) 



Im7r M y =tanh|- J 9g-tt|Ini7rn /lw 

_ »7ri2^e(g ■ u) 
2n~ j (x) 



(2.5) 



where we have defined 



Re 71"^ = -{tt^u + 71-*^) 



(2.6) 
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with similar definitions for Heirn^ and Imnn^. In Eq. ( |2.5| ) j3 is the inverse temperature, s(z) = 9{z) — 0{—z) and 
where, for the photon, 

x = /3q-u. (2.8) 

In our calculations we have adopted — (1,0), so that all the three-dimensional quantities refer to the frame in 
which the medium is at rest. 



p + q 




P 



FIG. 1. Diagram for the contribution to the photon self-energy matrix from a generic fermion. 

In practical calculations it is much simpler to determine the imaginary part of the self-energy using the second form 
in Eq. (pT5|), and therefore we will determine tt^P from the formulas 

Reir^P = ReTTn^ 
^ v 2n 1 {x) 



(2.9) 



The various components of the photon self-energy matrix are calculated from the diagram of Fig. [j]. For any fermion 
/ in the loop we use the propagator 



In these formulas 



with 



and 



^ ] 22 {p) = {i> + m f ) 



1 



jr — mi + ie 
-1 

p 2 — nij — ie 



+ 2TriS(p 2 — m 2 )r] / (p) 
+ 2iri5(p 2 ~ m 2 )r)f(p) 



Sf 12 (jp) = {£> + m f )2mS(p 2 - m 2 ) [ v ,(p) - 9(-p ■ u)\ 
S ( Jl(j>) = {£> + m f )2mS(p 2 - m 2 ) [ Vf (p) - 9(p ■ u)] . 



Vfip) = Q(P ■ u ) n F{Vf) + 0(-p ■ u)n F {~y s ) 



n-F(y) 



1 



e y + 1 



y f =Pp-U-a f , 



(2.10) 
(2.11) 

(2.12) 
(2.13) 



where a/ is the fermion chemical potential. 



3 



We c ons ider the calculation of ttu^, from which the real part of the physical self-energy is determined by means 
of Eq. (U). Referring to Fig. @, the contribution from one fermion in the loop is 



Jf) 



d4 P .•Hujo!/) 



•u^ = (-i)H) 2 Tr y ^ J i)7'(ff)isS l ( P + ?)iJr ) (-9WA(p) 

The electromagnetic couplings are defined by writing the matrix element 

.(err 



:(/) 



(f( P ')\jjf m Ho)\f(p)) = u(p')i ( ; u m) (9) u (f) - 



(2.14) 



(2.15) 



where q = p — p' is the momentum of the outgoing photon, u{jp) is a Dirac spinor and jfl^iq) is the total electro- 
magnetic current of the fermion. For the electron Jeps (?) = ^7^, while for the nucleons 



In Eq. (|t|) the anomalous part of the nucleon magnetic moments, given by 



(2.16) 



k p = 1.79 



K n = -1.91 



2m p 
( le 



\2m n 



(2.17) 

e stands for the electron charge and, as usual, a^u — ^ [7^,71,]. For arbitrary values of q, Eq. ( [2.16j ) does not hold 
and, in general, we have to parametrize the nucleon electromagnetic current in terms of the momentum-dependent 
form factors. However, in the applications of interest to us the photon momentum is sufficiently small that the form 
factors are well approximated by their values at q 2 = 0, and therefore Eq. (2.16) is appropriate. 
Substituting the formula for Sp^ given in Eq. ( 2.1C| ), fli{^„ can be written in the form 

(2.18) 

where iTii^ is the standard vacuum polarization term, while the background dependent contribution is given by 

d A p 



_(/) __.(o) , '(/) 



7V' (/) 



rTr 



3 ( f 7\q){£>+ h + m / )jjr ) (- 9 )(i» + m f ) 



(2tt)4 

ri f (p)6(p 2 -m 2 ) Vf(p + q)6((p + q) 2 -mj) 



(p + q)' 2 



rrij + it 



p 2 — m 2 + it 



2nir]f{p)r] f (p + q)5(p 2 - m 2 )S((p + q) 2 - mj) 



(2.19) 



The term in Eq. ( 2.19| ) that contains two factors of r/f contributes only to the imaginary part of the self-energy. We 
restrict our attention for the moment on the real part, which receives contributions only from the terms linear in r/f, 
and which we denote by Re T^t^v • Then, making the change of variables p + q — > p in the second of these terms in Eq. 
( p. 19 ) and carrying out the integral over poj we obtain from Eq. (2.G) 



R C7r (eff) _ _(7) + -(oO + -(no) 



(2.20) 



where 



T (7) 



-4e1 



d 3 P 



(27r)32£ (// + / / ) 
^ ) -- 4K //(2^^ + / 7 ) 



L, 



q 2 + 2p ■ q 



(27T) 

7r iX' j) = -8(e/K/m/)g 2 ff M!/ 



M„ 



q 2 + 2p ■ q 



(q 



(q 



-q) 



d 3 V 
(2n) 3 2£ 



(ff + fj) 



1 



q 2 + 2p ■ q 



+ (q^ -q) 



(2.21) 
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In these formulas, 



Lpv = 2p fl p L , + p^q v + q^pv - g^p ■ q , 

= q 2 (2m 2 - p ■ q)g^ - 2(p • qfg^v + 2(p • q)\PnQv + q»Pv] - 2q 2 PnPu , (2.22) 



P f 
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and ff j denote the particle and antiparticle number density distributions given by 

1 



S /3£T"/ + 1 



(2.23) 



(2.24) 



with the minus(plus) sign holding for the particles(antiparticles), respectively. The integrals in Eq. (2.21) is to be 
interpreted in the sense of its principal value part. 

As we have already mentioned, for the imaginary part of tt^u it is simpler to calculate Tti2^v first and then use 
the second form of the formula for Ixmrffl^ given in Eq. (2.9), instead of calculating Im7r llM1/ . Nevertheless, we have 



carried out the calculation in both ways and explicitly verified that they yield the same result for 7r^f . The results 



of that calculation are summarized in Section [V F , and the main steps leading to it are oulined in Appendix [B|. In 
the remainder of this paper we focus on the real part of it0 ■ 

In an isotropic medium, the most general form of the photon self-energy is [^|Jl3|| 



r (e#) = 



(2.25) 



where 



- Q 



v 



(2.26) 



with 

?V = g^u 1 - 

In general, ttt,l,p are functions of the scalar variables 

Q = q ■ u 



g = v^ 2 - q 2 



(2.27) 



(2.28) 



which have the interpretation of being the photon energy and momentum (in the rest frame of the medium) . The 
conditions under which the function irp can have a nonzero value, as well as some of its implications, have been 
discussed in Ref. In the present case, it is easily seen that it is zero since all the terms in Eq. (2.21) are 



symmetric in fi. v. The functions ttt,l are determined by projecting Eq. ( 2.2C ) with the tensors R^ and Q^. This 
procedure then leads to 



7T« = 

^ = 

fAV 



~4e 



u 2 
B' 

A' / 

/ ~2 
1 U z 



_ 2 Q MV 



-8(efKfm f )q 2 Df [R^ + Q^} 



(2.29) 



where, 
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Bf = 
Df = 



d 3 V 
{2tt) 3 2£ 

d 3 V 
(2tv) 3 2£ 

d 3 V 
(2tt) 3 2£ 

d 3 V 
{2tt) 3 2£ 

d 3 V 
{2tt) 3 2£ 



Us + fj) 
(// + fj) 
(// + fj) 
(// + //) 
(// + //) 





1 Or, . n 

■+■ zp q 




m u T 

U J-'/jtv 


a 2 


4- On ■ n 
■f if) q 




lvl n 


Q 


+ 2p-q 




u v M^ 


q 2 


+ 2p-q 




1 


q 2 


+ 2p-q 



+ (« 

+ (q 
+ (Q 
+ (Q 



-q) 
-q) 
-q) 
-q) 
-q) 



(2.30) 



Substituting here the explicit formulas for and given in Eqs. J2.22| ) and ( |2.22| ), after some straightforward 
algebra we obtain the relations 



A' 



T 1 \a } + 3mjD f 
1 



B' f = -Q 2 [-A f + m 2 D f \-q 2 B f 



(2.31) 



with Af,Bf being given by 



A, 



d 3 V 



B f = 



(2tt) 3 2£ 

d 3 V 
{2n) 3 2£ 



(// + fj) 
iff + //) 



2m 2 : — 2p ■ q 
q 2 + 2p ■ q 



+ {q^ -q) 



2(p ■ v) 2 + 2{p ■ v) (q ■ v) - p ■ q 
q 2 + 2p ■ q 



+ {q^ -q) 



(2.32) 



For l ater p urposes, it is also useful to note here that, by inspection of the formulas for Df and Af given in Eqs. ( [2~30| ) 
and (2.32) respectively it follows immediately that 



where 



A f = -2J f + (2m 2 f + q 2 )D f 
Jf = J (2^ (// + / 7 } - 



(2.33) 



(2.34) 



From Eq. (2.29) we can immediately identify the contribution of any fermion to the transverse and longitudinal part 

( f) A 1\ 

of the self-energy, which we denote by tt^ l . Using Eq. (2.31) they can be expressed in the form 



Re 7TjP = -2e 2 f ( .1, + -^IJ 

Re Jf)^A„ 2 l 



Q' 2 



8efKfmfq 2 Df — 2n 2 q 2 



2m Pf ~ ^ B f 



4<Sf Q2 B/ ~ 8e/K/TO/<7 Df ~ 4K f q 



2 A f + ^ B f+ m f D f 



(2.35) 



In this way, the contribution of any fermion to the photon self-energy is expressed in terms of the three functions 
Af,Bf,Df. While the evaluation of these functions is not possible in the general case, some useful results can be 
obtained by considering special cases. Of particular interest to us are the approximate expressions obtained for small 



values of q, which we analyze in some detail in Section [II 



III. LOW PHOTON-MOMENTUM LIMIT 



The functions Af,Bf,Df are the same ones that appear in the calculation of the induced electromagnetic vertex of 
a neutrino in a matter background, which were analyzed in Refs. |l4| , [i1| with considerable detail for various limiting 
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values of the photon momentum, and for several conditions of a background gas made of nucleons and electrons.^] 
The results obtained there are directly applicable here. In this section we expand upon those results further. 

In most of the situations of practical interest, the photon momentum is such that CI, Q <C m p ^ n . In this case, from 
Eq. fl2.33|) , it follows that 



D, 



2m 2 f 



(Af + 2J f ) 



(3.1) 



Later on we use this equation to evaluate D n ^ p . On the other hand, borrowing the results given in Eqs. (A6) and 
(All) of Ref. (l5|, for small values of q we have 



B, 



A, 



d 3 V 



v v ■ Q 



B f + - 



(2tt) 3 \n-v-p-QJ d£ 



n f d 3 v 



(2n) 



(ff + fj)- 



(3.2) 



with v-p = V jCl denoting the velocity of the particles in the background. The above formula for Bf can be rewritten 
by multiplying the integrand by the factor 



-(Cl-v-p-Q + v-p-Q). 
The first two terms integrate to zero, while the third one leads to 



B, 



1 



d 3 V ( {v v -Qf 



(2tt) 3 \ n-vvQ d£ 



(If + fj) 



(3.3) 



(3.4) 



We should note that expressions in Eq. (3.2) are derived from Eq. ( 2.32j ) by expanding the integrands in terms of q/£ 
and retaining only the terms that are dominant when q/£ — > 0. Therefore, Eq. (3.2) is valid for values of q such that 



q/(£) « 1 



(3.5) 



where (£) denotes a typical average energy of the fcrmio ns i n the gas. Thus, for a non-relativistic gas, Eq. (3.2) holds 
for f2, Q <C to/. If the gas is extremely relativistic, Eq. ( |3.2[ ) also holds for fl, Q > mf, subject to Eq. fl3.5|). 
For an isotropic fermion distribution, an elementary integration on the angular variables yields 



B, 



A, 



1 

Bf + 



dVV*- {ff 



h 



£n 

2TQ 



In 



£n + VQ 



4tt 2 



dV 



r^d_ 

~~£d£ 



(// + fj) 



2Q 



In 



£n-pQ 

£n + VQ 



£il-VQ 



(3.6) 



Up to this point no assumption has been made regarding the nature of the fermion background. Apart from the 
restriction on q, Eqs. (3.2) and (3.6) hold for a relativistic or non-relativistic gas, whether it is degenerate or not. 
Accordingly, Eq. (3.6) serves as a convenient starting point to find the dispersion relations of photons that propagate 
through an isotropic medium, provided Eq. ( |3.5| ) is verified. In what follows, we present the explicit results for several 
important physical situations. 



A. Degenerate gas 

As functions of the energy, the distribution for a degenerate gas and its derivative are given by 

ft = 9{£ F f - £) , 



(3.7) 



1 The quantities denoted by A e , B e in the present paper were denoted by A,B in Refs. 
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As a consequence, the integrations needed to determine Af,Bf,J~f become very simple. The results can be expressed 
in the form 



Jf = 
Bf = 
Af = 



ra s V Fflf 
8tt 2 

4tt 2 

m f V Ff jf 
Air 2 



± + VFf 

2j'jv F f "° V 1 - l 'Ff 



■log 



-1 + 2 Z f lo S 



l + z f 



1 + Zf 



1-*/ 



(3.8) 



where 



7/ 



l-v Ff 



VFfQ 



with upy = Vpf/^Ff being the Fermi velocity of the gas. Let us note that 

A f (Q=Q) = -2J f , 



(3.9) 



(3.10) 



which is va lid in gener al, fo r any distribution function and any kinematic regime of the gas, as can be easily checked 
from Eqs. (|3|) and (OJ). 



B. Relativistic gas 



If the fermions are ultrarelativistic, then £ ~ V and Eq. (3.6) reduces to 

n + Q 



U f = I 1 - — In 



V 

Af = —3u>, 



2Q 



of ■ 



The quantity u)Qf is given by 



-o/ 



d 3 P 



{2nf2£ 



(ff + fj) 



1 - 



Z£ 2 



and in Eq. (3.11) we have used its limiting form 



,J of 



^f d VV(f f + f 7 ), 



for £ = V . In terms of Eq. ( 3.13|) Jf becomes 



2 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



The remaining integral in Eq. ( |3.13| ) cannot be performed without specifying the distribution function. For a 
classical (relativistic) gas // = exp(—(3fV + aj) and 



2 Pf 



while for a degenerate gas the integration over V yields 



1 /3nA 2/3 = V\ t 
12 V 7T / 12tt 2 



(3.15) 



(3.16) 



The formulas for Af and Bf given in Eq. ( 3.11 ), with determined by Eq. ( 3.16 ), coincide with the relativistic 
limit (£p — Vf) of Eq. (3.8), as it should be. 
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C. Non-relativistic, non-degenerate gas 



The integrals cannot be carried out explicitly in this case. However, some useful approximate results can be obtained 
by considering the two regions O 3> zJ/Q and fl <^TJfQ separately, where 



v 2 ~ 



1 



/3 fm f 



(3.17) 



We have denoted the inverse temperature of the gas by f3f in order to allow for the possibility that different components 
of the background may be at different temperatures. 

The results that we have obtained in this way can be written in the form 



Bf = 
Af = 



f n f J ~} + z2 f ( for z f < 1) 



4 

Pf n f 
4 

"/ 

4m f ' 



+ A (for z f > 1) , 

L . 3 —2 



(for Zf < 1) 
+ (forz / »l). 



with Vf as defined in Eq. ( [3.17 ), and 







(3.18) 

(3.19) 
(3.20) 

(3.21) 



As an illustration of how we have proceeded, let us consider in some detail the calculation of Bf. For a classical 
non-relativistic gas, 



and from Eq. (3.2) 



where 



B I 



d£ 



(3 f n f n(3 f m f 



-Ms » 



dVVff log 



vvQ 



For the low frequency region, we proceed by putting 



log 



and then using the result 



dVff 



(3.22) 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



in Eq. ( [3.23] ). In this way the result for <C VfQ given in Eq. ( [3.18 ) follows. For the high frequency region, we put 

1 1 



log 



(3.27) 



and then use, in addition to Eq. ( 3.26| ), the results 
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/ dVV 2 f f =7T 2 n f , 
Jo 

Jo PS 

/ dVP 6 f f 
Jo 



157r miuf 
ft 



in Eq. ( p3| ). 

Similarly for Af, following a similar procedure it is easy to show from Eq. (3.2) that 



n S 



4m 



(^) 2 (forfi»T7 / Q) 



/ 3 + 5 



Therefore, comparing with the formula for Bf in Eq. (3.18), we arrive at Eq. ( |3~T9| ). 



(3.28) 



(3.29) 



IV. DISCUSSION AND APPLICATIONS 



In ge nera l, the dispersion relations are obtained by solving Eq. (2J) with — which, using the decomposition 
in Eq. ( [2.25| ) and the fact that up = in the present context, is equivalent to solve 



[iV(-<Z 2 + ttt) + Q^u{-q 2 + 7T £ )] A v = . 



(4.1) 



Since R^ v and Q M „ are orthogonal, the solutions to this equation exist only for q = (uij \l{Q) , Q), where lot,l(Q) 
satisfy 



_ Q 2 = tt(ut, Q) i 
u 2 L -Q 2 =tt l (lo l ,Q). 



(4.2) 



In general, lot,l are complex functions of Q and can be written in the form ujt,l = ^t,l — «7t,l/2. The quantities CIt,l 
and 7t,l are real and have the interpretation of being the dispersion relation and damping rate of the propagating 
mode, respectively. Retaining terms that are at most linear in 7t,l it follows that 



T,L 



Q 2 =Rei: T . L (n TiL ,Q), 



1T,L 



Im ttj 1 j 



R.e7r Tji ) 



(4.3) 



For a background containing several fcrmions species, 



(4.4) 



Each individual contribution in Eq. (4.4) can be immediately identified from Eq. (2.29). Thus, for the electron 
component we set n e = while for the neutron e n = 0. The corresponding expressions for the proton are more 
complicated due to the presence of two terms in jf™^ ■ 

The explicit formulas given in Section III allow us to study in some detail the photon self-energy in a nuclear 
medium, provided of course that the physical conditions are such that the approximations and idealizations that led 
to them are valid in the particular context in which they are being applied. As an illustration, we consider some 
possible situations. 



10 



A. Electron background 



Although the results for this case ar e we ll known, we briefly review them here in a form that will be useful for our 
later purposes. Setting n e = in Eq. ( 2.35 ), 



Re4 6) = 



-2e 2 A 



Re 7r 



(e) 



4e 2 — B 



(4.5) 



The functions A e and B e are in given in general by Eq. ( 2.32 ) and, in the low q regime, by Eq. (3.2). By specializing 
these formulas further to the zero frequency (static) limit or zero momentum (long wavelength) limit, we obtain explicit 
formulas for them, given explicitly in Refs. jl4 15 1, which are useful in practical applications. For example, from Eq. 
(3.9) of Ref. g5j we have 

^(fi,Q-»0)=-3wg e> 

,2 



a 2 



(4.6) 



Substituting Eq. (4.6) into Eq. (4.5) we then obtain 



Re 7r, 



, T)i (f2,0)=4e^ e , (4.7) 

which in turns implies the well known result that the longitudinal and transverse permitivity have the same long 
wavelength limit 



ew(n,o) = 1 



4^4 

n 2 



(4.8) 



Eq. (4/7), or equivalently | 
namely 



implies the well known expression for the electron contribution to the plasma frequency, 



Ae 2 u 



0. 



defined as the value Ot,l(0) of the dispersion relation in the limit Q 

nonrelativistic, then uj 2 e — n e /4m e , as is imm ediately seen from Eq. ( 3.12 ), and fl e — \J e 2 n e / 
regime, u> 2 e is given by Eqs. ( 3.15| ) and ( |3.16|). 

For values of Q ^ 0, using Eqs. (U) and (j3j) in Eq. @ we obtain 

d 3 r 



(4.9) 

In particular, if the electrons are 
In the relativistic 



Re7ri 



Re7r ej =- 



e 2 n 

2eV 

Tic? 



(2tt) 3 
d 3 P 



where 



(27T) 



V-p 



V± ■ V-p 

( (vv ■ Q? ^ 



{vv ■ Q)Q. 



d£ 



(fe + fe) ■ 



dg (fe + fe), 



(4.10) 



(4.11) 



Substituting Eq. (4.12) into Eq. (A4) we obtain the explicit formulas for the transverse and longitudinal components 
tt,i of the dielectric constant, that coincide precisely with those obtained by the semiclassical approach based on the 
Boltzman equation Jig ], For the particular case of an ultrarelativistic electron gas, using the result in Eq. ( 3.1l| ) we 
obtain 

„2 



IW/' = * I ^-^ln 

2Q 



fl (2VL 



n + Q 



n-Q 
n + Q 



(4.12) 



Another useful limit corresponds to set m e — > in Eq. ( [2.32| ). The formulas obtained for Re7r^ L by substituting 



the expressions so obtained for A e and B e into Eq. ( |4.5| ) , coincide with those given by Weldon [|2J . Those formulas are 
valid in the regime T,q 3> m e , and neglecting the terms of order q/T they coincide precisely with the result given in 
Eq. ( 4.12|) . In contrast, the formulas given in Eq. (^2) and the corresponding expressions for Re7r^ L obtained from 
Eq. (4jf) are valid for q <C {£), but they hold whether T is larger or smaller than m e . 
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B. Degenerate neutron gas 



The neutron con tribution to the photon self-energy is given by setting e n — in Eq. ( 2.35| ). Remembering Eq. (3.1) 
and using Eq. (|3.8|), we obtain for this case 



T) ( n ) 2 M 
Reir y T — q Xq In 



1 

2^ 



1 



Re 7r 



(n) 



2 (n) 2 
Q Xq lnV Fn 



2v F 

1 



log 



Kn 



1 



1 + VFr, 
1 - «Fr, 

log 



- 1 



1 + VFn 
1 - V Fn 



where we have defined 



1- 2 W f„(1 + 4) 



1 - -z„ log 



+ 4+ 2 Z «( 1_Z n) l0 e 



(„) K^mnVFn 
Xo = "a 



1 + z„ 



1 - z„ 



1 + z„ 



(4.13) 



(4.14) 



We emphasize that the above formulas are valid for all values of q subject only to fi, Q <C m n , and also they hold in 
the non-relativistic as well as the relativistic case. However, in the non-relativistic limit, they reduce to 



Revr^ = q 2 Xo n) fT(z n ) , 
Re7r[ n) = q 2 Xa n) v 2 Fn fL(zn) , 



(4.15) 



where 



-I ± 2 2 
1 - -V W Z 



1 z log 

2 B 



1 + z 



i + z 2 + iz(l-z 2 )log 



1 + z 



1 



(4.16) 



In the limits — ► or Q — > we obtain from Eq. (4.15) the limiting values 



Re4")(0,Q) = -Q^, 

Re4 n) (o,e) = ie a ^ n) <, 



R e 4 n) (n,0) = Re7rr(fi,0) = ^X«„- 



» 



l2,i«)„,2 



(4.17) 



I n pa rticular, it is instructive to observe that taking the static limit (Q = 0) in Eq. (A8) and using the results of Eq. 
([4.17D , we obtain 



M(o,e) 



(n) 

X 



(4.18) 



which, glanzing at Eq. ( 4. 14] ) , is recognized as the classic Pauli formula for the static magnetic permeability of a 
degenerate fermion gas. For later reference, it is useful to observe that in the nonrelativistic limit 



x[, n) ~ (3.6)av Fr , 



(4.19) 



which follows from Eq. (4.14) by using Eq. ( 2.17] ) and setting VFn ~ fn n VFn- 

For our purposes, the primary quantities of interest are ttt.Li which determine the dispersion relations of the photon 
mode s through Eq. (4.2). If the electron and proton contribution to the self-energy can be neglected and we use Eq. 
( 4.15 ) for the neutron contribution, the dispersion relations are obtained by solving the equations 



(n 



1 - Xo JT 



tin 



rn? - q 2 



1 Xo v FnJL 



VFnQ 
VFnQ 



= 0. 
= 0. 



(4.20) 



The functions /t,l are ploted in Fig. ||. 
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"1 1 — 

f T ,v F = 0.1 
f T ,v F = 0.3 
h 



2 4 6 8 10 

z 

FIG. 2. Plot of the functions ft and fx defined in Eq. ( [4. 16) ) . At z = 1, fr becomes infinitely negative. Even for relatively 
large values of vf the function fr has approximately the same form, and for small values of vf the curves fall almost on top 
of each other. 



As can be seen from the figure, both functions satisfy Jt,l < 1- On the other hand, from Eq. (4.19) we have 

Xo < 1. Therefore Eq. ( 4.2C ) has no other solutions than Qt.l = Q, so that the photon dispersion relation is not 
modified by the presence of the background in this case. 



C. Relativistic electron and degenerate neutron gas 



In general, the total background contribution to the transverse and longitudinal components of the photon self- 
energy in a matter background is 



0) , (p) 



In) 



(4.21) 



Let us neglect for the moment the contribution from the protons, and assume that the electrons c an be represented 
by a relativistic gas and the neutrons by a degenerate non-relativistic gas as considered in Section IV B . In this case 
7Ty ^ are given by Eqs. (4.12) and (4.15), respectively. 



The dispersion relations for the transverse modes in this case are obtained by solving the equations 



g 2 =Re4 n) 



RC7T 



(e) 
T > 



(4.22) 



which can be written in the form 



<1~ 



1 - Xo JT 



o 2 2 *" 
= 3e W 0e^ 



2tt q 2 , 



fi + Q 



(4.23) 



While it is not possible to find the general solution to this equation, some useful conclusions can be drawn from it. 
For example, for values of < Q, the right hand side of Eq. (4.23) is a positive quantity, while the left hand side 

Therefore, the solution to Eq. (4.23) is such that fix > Q- 



is negative since, as we have seen above. 



Xo 



It < 1. 



This implies, in particular, that a particle propagating such a medium cannot emit transverse photons in the form 
of Cerenkov radiation. Furthermore, since z n is considera bly g reater than unity for CI > Q, the value of /t(zti) is 
negligible for that range and therefore the solution to Eq. ( 4.23| ) is well approximated by the corresponding solution 
for the electron background only. 

Notice that a possibly different, and erroneous, conclusion would have been obtained if we were to approximate the 
neut ron c ontribution in Eq. ( 4.22 ) by the value in the static limit ReTr^^O, Q) given in Eq. ( 4.17 ). Then, instead of 



Eq. (4.23), such an approach leads to the equation 



-Q 2 x { o 



Re7r 



(e) 



(4.24) 



for the dispersion relation, which can have a solution such that Qt < Q depending on the relative size of the two 
competing terms in the right hand side. 
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For the longitudinal modes on the other hand, the dispersion relation is determined by solving the equation 



with Re?r[ e) given by Eq. ( [4.121) . The factor 



n 



n + Q 



n-Q 



(4.25) 



(4.26) 



which is plotted in Fig. ||, is larger than unity for values of Q/Q larger than about 0.84. Therefore, there is a small 
range 



< n 

0.84 ~ — < 1 



Q 



for which the solution to Eq. (4.25) is such that Ql < Q 




0.2 0.4 0.6 0.8 1 

n 
Q 

2Qi««r|n±fi 



FIG. 3. The factor ^ log | f|±j| | is plotted as a function of ^. 



(4.27) 



Hence, the Cerenkov radiation of longitudinal photons is allowed. For those values of f2/Q, the value of z n corre- 
sponds to the region in Fig. || where the function _/x is negligible. Therefore, the neutron component of the background 
does not play a significant role in the photon dispersion relation in this case either. 



D. Relativistic electron and degenerate neutron and proton gases 

The contribution from the protons can be obtained straightforwardly by substituting Eqs. ( |3.l| ) and ( ft.8| ) in Eq. 
(2.35). In the non-relativistic limit, the result can be expressed as 



Re4 p) = 6eW 0p g T (z p ) + q\^ f T (z p ) + q 2 (-^-) X ^h(z p ) , 

\ Trip hip J 

Re^ = UeW Qp9L (z p ) + q\^ f L (z p ) + q 2 (-^-) X ^h(z p ) , 



where, in analogy with Eq. (4.14) 



Xo 



( p ) _ n 2 m p V Fp 



and 



0p 4m p ' 



(4.28) 



(4.29) 



(4.30) 



with n p denoting the proton number density. The functions Jt.l are defined in Eq. (4.16), while 
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g T {z) = z 2 + -z(l-z 2 )log 



1 + z 



9l(z) -- 
h{z) = 

These functions are plotted in Fig. |]. 

2 
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1 



1 -z 
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(4.31) 



FIG. 4. Plot of the functions h and gp,L for some representative values of vp. At z = 1, <?l,/i become infinitely negative. 
For 2> 1, the functions h,gp,gL have the asymptotic values ~u^,2/3, 1/3, respectively. 



Let us consider the dispersion relation for the transverse modes. It can be written in the form 

q 2 (l - x) = Re4 e) + 6e 2 io 2 p g T (z p ) 

where we have defined 



(p) 



» 



Trip fop 



h(z p ) . 



(4.32) 



(4.33) 



For £1 < Q, we have tt^ > 0, as we al so ar gued before in Section [V C , and from Fig. [| we can see that gp > also. 
Therefore, the right-hand side of Eq. (4.32) is a positive quantity. On the other hand, since the functions h and fr 



never become l arger than unity, then \ < 1 since Xo™'^ are considerably smaller than unity (they are of order avp). 



Therefore, Eq. ( J4.32) ) has no solution with Q < Q. For the opposite regime > <2, we have z UiV ^> 1, and fp,h,gT 
go to their asymptotic values 



9t 



J F . 



J F . 



(4.34) 



It then follows that the contribution from the \ terms in Eq. (4.32) is small. On the right-hand side of that equation 



the term Wq p 5t from the protons is competing against the corresponding term from the electrons. Since 



J 0e 



and n e ~ T 3 , where T is the temperature of the electron gas, we have oj 2 e 
from Eq. ( |43C| ) it follows that 



n p /T, assuming that 



T 



J 0p 



n p . Then 



(4.35) 



so that the proton term on the right-hand side of Eq. (4.32) is also small. The equation for the dispersion relation is 
then well approximated by simply 



2 (e) 

q =ir y T > 



(4.36) 
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Thus, the nucleon contribution, represented by the x an d <?t terms in Eq. (4.32), may affect some of the details of 
the transverse dispersion relation for the relativistic electron gas, but it will not modify the general features in an 
essential way. 

The situation can be different for the longitudinal modes. The equation for the dispersion relation can be written 
in the form 



q 2 (l - x') = Re^ e) + 12e 2 uj 2 p g L (z p ) 



(4.37) 



where 



I — (n) .2 
X = Xo V F; 



fdz n )+x^v 2 Fp fL(z p ) + X^ 



(P) 



h(z p ) 



(4.38) 



The factor v 2 F in the terms proportional to /z,(-Zn,p) makes those terms unimportant, while the term proportional to 
h nee d not be negligible, although it is small as discussed above. Regarding the terms on the right-hand side of Eq. 
(4.37), consider first the region > Q. In that region, we have z p ^> 1 and hence gL — > 1/3. The same argument as 
before [Eq. (4.35)] then implies that the proton contribution is not relevant. However, in contrast with the equation 
for the transverse modes, Eq. (4.37) can have solutions with Ql < Q because <?z, can be negative and, as we saw in 



Section IV C, tt^ can be negative as well in this range. Furthermore, for values of fi, Q such that z p ~ 1, we have 
ql ~ — -2-, and this can bring the proton term uJn v gL to be comparable to the electron term wL Thus, to summarize 

V F J- 

the situation for the longitudinal modes, in Eq. (4.37) the term with fa is irrelevant, the term with h can give a small 
correction, while the term with g^ can have an important effect in the dispersion relation. The latter effect can occur 
in the range of Q, Q such that z p ~ 1, if 



,2 



■ J 0e 



(4.39) 



Using Eq. ( 4.35 ) this can be translated into the condition v Fp ~ where T stands for the electron gas temperature. 



E. Relativistic electron and classical neutron and proton gases 



For a classic al n on-rel ativis ti c pro ton gas, th e con tribution to the self-energy is obtained by inserting the expressions 
given in Eqs. (3T) and ( 3.18 - 3.20| ) into Eq. ( 2.35| ). For the transverse part we thus obtain 



Re 4?) = \A^L Gt {z p ) + 4n p P p q 2 F T (z p ) + \^n p (3 p q 2 H(z p ) 



where, to the order that we have calculated, 

G T (z f ) 

F T = H(z f ) 

Similarly, for the longitudinal part, 



z) (z f « 1) 
1 + ^ (2/>l) 

1-zj (z/«l) 
TT 2 -i (z/»l) ' 



Re7r ( p)= \efn P 



G L (z p ) + ^ q 2 F L (z P ) + ^n p P p q 2 H{z p ) . 



where 



G L (z f ) 
F L (z f ) 



^\z 2 ~ 



-1 + zj (zf <C 1) 



-l + 2z^ (zf < 1) 
1 + ^ (z f » 1) 



(4.40) 



(4.41) 



(4.42) 



(4.43) 



1G 



The formulas for the neutron contribution correspond to setting |e| = in Eqs. ( [1.40 ) and (1.42) and replacing every 
quantity that refers to the proton by the analogous quan tity for the neutron. 

Comparing these results with those given in Eq. ( 4.2§| ) for the degenerate case, it is easily seen that the qualitative 
features betwee n that cas e and the present one are similar. In fact, an analysis along the same lines to that carried 
out in Sections IVC and IV Ej reveals similar conclusions. In particular, the neutron component does not modify 
the photon dispersion relations in a significant way. The proton contribution to the transverse component of the 
self-energy is also unimportant but it can have a noticeable effect for the longitudinal part. This can occur for values 
of il < Q, but such that z p ~ 1. In this case, similarly to the situation for the case of degenerate protons, Gl ~ l/v p , 



and this brings the proton contribution to Kenr, to be of the order of e n p [3 p . Assuming that n p 



and rip ~ T 3 , 



where T stands for the electron gas temperature, the proton contribution is then of order e 2 T 3 /3 p . This can compete 
with the electron contribution, which is of order e 2 T 2 , if the temperature of both components are comparable. The 
proton contribution could in fact dominate if the proton gas temperature is smaller than that of the electron gas. 



F. Nucleon contribution to the imaginary part 



In general, the concept of a propagating mode with a definite dispers ion relation is meaningful provided that the 
damping rate 7t.l is much smaller than £It,l in the solution to Eq. (1.2). Therefore, in the cases in which the 
nucleon contribution to the real part of the photon self-energy may be substantial, it is imperative to know also the 
contribution to the imaginary part. 

In the general case, the calculation of the nucleon contribution to the imaginar y pa rt of t he s elf-energy is cumber- 
some. Fortunately, for our purposes, in which the small q formulas given in Eqs. ( |3.l| ) and (3.2) hold, we can give a 
simple recipe to determine it. It is given by the so-called Landau rule, whose justification in the present context is 
provided i n A ppendix 

from Eq. (^2) by making the replacement f2 
part of the resulting expressions. Thus, 



The rule is simply that, in the small q limit, the imaginary part of Af,Bf can be obtained 

+ it in the integrand of those formulas, and taking the imaginary 



d 



ImAi 



lmD f = 



Further, the imaginary part of t^tl ^ s gi yen by formulas analagous to Eq. (2.35), but with the functions Af,Bf,Df 



2m 2 f 



ImAf 



.(l-v 2 v )5(Q-vvQ) — (ff + fj. 



(4.44) 



replaced by their imaginary part, given in Eq. (4.44) 



1. Degenerate case 

For the particular case of the degenerate gas, these formulas reduce to 



ImB 



/ 



ImAf = 2m 2 f lmD f = (1 - Vp f )lmB f 



z f 0(l-\z f \), 

Ff) 



where 7/ and Zf are defined in Eq. (|3.9|). These in turn yield, 



(4.45) 



Im 7r 



(/) 



" \zf\hf < Se 2 tu 2 0p (z 2 -l) + q 



Jmni f) =nz f 6(l-\z f \h A -6. 



2m,fK p j 2 




q 2 X ^v 2 F (z 2 -l) 



(4.46) 
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2. Classical, non-relativistic case 



The formu las in Eq. (4.44) can be evaluated explicitly also for the classical non-relativistic case. With z; as defined 
m Eq. the results can be expressed in the form 



ImB 



1 hz 



f 



4 V 2 



P f n f z f e ? z ! 



ImAf = 2mj lmD f = ImB/fl - vj(2 + zj)] 



(4.47) 



which imply 



Im 7r, 



(/) 



-(3 f n f z f e-^f i -ejv} + ^% 2 (1 - vjzj) + k)<?{\ - vjzj) \ , 



Im 7r 



(/) 



f-P f n f z f e-^ {e}(l-v}z}) + ^(1-vjzj) - 2^} . 



(4.48) 



As usual, these formulas allow us to obtain the results for the proton and neutron separately, by setting e / and Kf to 
the appropriate value in each case. 

V. CONCLUSIONS 

In this work we have calculated the 1-loop contribution to the photon self-energy in a background composed 
of electrons and nucleons, taking into account the anomalous magnetic moment couplings of the nucleons to the 
photon. The generic contribution from a fermion to the real part of the photon self-energy is given by the formulas 
in Eq. (2.35), in terms of the quantities that we have denoted by Af,Bf,Df. The explicit evaluation of these three 
functions is carried out in Section [II for various limiting cases of interest. The results obtained there were applied 
in Section IV to consider the effect of including the nucleoli contribution to the photon dispersion relation in a 
plasma. For concreteness, we considered the case of a relativistic electron gas, superimposed on a nonrelativistic 
nucleon background, which we considered to be degenerate or classical in separate cases. The results of that exercise 
indicate that the magnetic moment couplings of the nucleons do not have a significant effect on the photon dispersion 
relations. Hence the neutrons, which only have these type of coupling, do not seem to play a siginificant role in the 
determination of the photon dispersion relations in such media. For the protons, the normal couplings do not have a 
significant effect on the transverse modes eit her, b ut t hey ca n have an impact on the longitudinal dispersion relation 
in some circumstances. As shown in Sections [VD and IV E, for values of the photon momentum such that Q ^v p Q, 



where v p is a typical velocity of the protons in the gas, the proton contribution to the longitudinal component of the 
photon self-energy can be of the same order as that of the electrons. Of course we have reached these conclusions by 
considering in detail the application of the general formulas to some particular idealized situations. In any case, the 
formulas that we have provided in this work open the way to similar studies in other possible settings. 

This work has been partially supported by the U.S. National Science Foundation Grant PHY-9600924 (JFN) and 
by Grant No. DGAPA-IN100694 at the Universidad Nacional Autonoma de Mexico (JCD). 

APPENDIX A: DIELECTRIC FUNCTION 

Although the properties of photons propagating through the medium can be deduced from the knowledge of the 
photon self-energy directly, it is useful to translate those results into the language of macroscopic electrodynamics. 
The transverse and longitudinal components of the dielectric constant of the medium are introduced by writing the 
induced current, in the rest frame of the medium, as 



find) = tn r (1 _ + (1 _ et) ^ + i£p Q x g 



where the longitudinal and transverse components of the electric field are defined by 



(Al) 



Ei = Q{Q ■ E) , 
E t = E — Ei . 



(A2) 
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This is the most general form of the induced current, involving terms that are linear in the field, and subject only to 
the assu mption of isotropy. On the other hand, by comparing Maxwell's equations with the effective field equation 
Eq. (2.1), it follows that the induced current is given by 



which can be written as in Eq. (|Al|) with 



Aind) 



l-e t 
1-ej 



7I"t 

9 ' 

T 

7Tp 

n 2 ' 



(A3) 



(A4) 



Alternatively to ey, the dielectric and magnetic permeability functions e,/x are introduced by writing the induced 
current in the equivalent form 



(1 - e)E + ill - -)qx B + ie, p Q x E 



(A5) 



where 



(A6) 



Comparing Eqs. (Al) and (A5), the relations 



e = ( 4 

1 i n \ 
M Or 



(A7) 



then follow, and in particular 



1 

- = 1 



fl 2 

ITT 5-7TL 

<7 Z 



(A8) 



APPENDIX B: IMAGINARY PART 



From Eq. (|2.3D , it follows that 



(Bl) 



where Im 7r M „ can be determined in terms of the components of the self-energy matrix by either one of the formulas 



given in Eq. (2.5). However, as already mentioned in Section ||, it is much easier in practice to calculate 7r 12Mi , than 
it is to calculate Im7rn M i/, although the procedure in both cases is similar. As an illustration, we outline below the 
main steps of the calculation in terms of i:i2^v , although we have verified explicitly that the same result is obtained 
using the alternate form in terms of Ivann^u. We start by considering the contribution from the electron loop, and 
then generalize the result to the case of the nucleons afterwards. 
From Fig. [j], 



(e) 



(-l)(-ie)(ie)Tr 



d 4 p 



(B2) 



where 



P =p + q- 

After substituting the propagators given in Eq. ( 2.10| ) we use the relation 



(B3) 



19 



faeO') - 0(-p' ■ v)][r) e (p) - 9(p ■ v)] = -n(x)[rj e (p)e(p' ■ v) - r) e (p')e(p ■ v) 

+ 6{p ■ v)9(-p' ■ v) - 9{-p ■ v)6(p' ■ v)] , 



which is proven by using the identities 



together with 



(B4) 

(B5) 
(B6) 
(B7) 

while x and y e are defined in Eqs. (2.8) and ( 2.1 3| ) , respectively, and we have used the fact that y' e = y e + x. In this 
way we obtain 



In Eq. (JB6J) we have defined 



rieip) - 9(p ■ v) = -e v "n F (y e )e(p ■ v) 
Ve(p) - 0(~P ■ v) = n F (y e )e(p ■ v) , 



e y 'n F {y e )n F {y' e ) = n 7 (x)(n F (y e ) - n F (y' e )) . 



y' e =p' -a e , 



™V2^ = 167r 2 e 2 n 7 (a;) 



d 4 p 



L^S(p -m e )5(p -m e ) 



x {r] e (p)e(p' ■ v) - T]e{p')e(p ■ v) + 9(p ■ v)9(-p' ■ v) - 9(-p ■ v)9{p' ■ v)} . (B8) 
Making the change of variable p + q — > p in the second and fourth terms in braces, Eq. (|2.9[) then yields 



Im n^ffj = ine 2 



d 4 p 



{L,J[{p + q) 2 -m 2 e )]5{p 2 -m 2 e ) 



x [ Ve (p)e((p + q)-v)+ 9{p ■ v)9(-(p + q) ■ v)] + (q^ -q)} , (B9) 
for the electron contribution to the imaginary part of tt^\ After doing the integral over p$, this reduces to 

d 3 P 



Im ir^J = Aire 2 



{L^5[{p + q) 2 -m 2 e 



(27r) 3 2£ 

x [9{{p + q) ■ v)(f e + /-) + 9(-(p + q). V )(l-f e - /-)] - (q - -q)} 



(BIO) 



Decomposing tx^J according to Eq. ( |2.3| ), Eq. ( B10| ) then implies 



where 



Im^L = 7r 



d 3 V 



Im7r ( r e) = 4e 2 ^ImB e , 
Q 2 



(2m 2 e - 2p ■ q)S[{p + q) 2 - m 2 e 



(Bll) 



Im Bp — 7r 



(2tt) 3 2£ 

x [9((p + q) ■ v)(f e + /-) + 9{-{jp + q) ■ - f e - /-)] - (q -g) 
d 3 P 



[2(p • v) 2 + 2{p-v)(q-v)-p- q}6[(p + qf - m{ 



(2tt) 3 2£: 

X [9{(p + q) ■ v)Ue + k) + 9{-{p + q) • w)(l - f e - h)] - (? - -q) 



(B12) 



The formulas in Eq. ( B12j ) compl eme nt those in Eq. (2.32), which hold for the real part. The counterpart of the 
approximate formulas given in Eq. (|3.2[) for the small q regime are obtained as follows. For q <C (£ ), we can use 



9[~(p + q)-v}S[(p + q) 2 -m 2 e }=0 



(B13) 
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to discard some terms in Eq. ( B12| ). The terms discarded in this way are the standard vacuum terms as well as the 
temperature corrections to them, which they are relevant only for q > 2m e and physically correspond to the creation 
of electron-positron pairs. The surviving terms can be written in the form 



ImB e 
Im A e 



d 3 V 



(2tt) 3 



{2£l + £sSl + V-q) 

-6{n + £p- Ep^) (f e + h) -(q^ -q) 



d 3 V J (2P 2 - UpQ, + V ■ q) 



(2tt) e 



AF -F, - 



5(n + £p-S p+? )(f e + h)-(q 



-q) 



(B14) 



where we have made explicit the dependence of the energy on the momentum variable. Making the change of variable 
V — ► V — 7:q and expanding the integrand up to terms that are linear in q, we obtain 



IrnB, 



-n 



d z V - d 

—W-V.q) — (f e + fe), 



(2tt) 3 



ImA e =ImB e + -SI 



d 3 V 
(2tt) 3 



v 2 v S(n-V-q)-^(f e : h 



(B15) 



The results given in Eq. ( B15 ), together with Eq. ( Bll ), provide the justification of the rule stated in Section 



■ (e)' 

IV F , as can be immediately confirmed by applying it to determine the imaginary part of it T L from the formulas 
in Eq. (4.1C) for the real part of the same quantities. The same arguments given here also to apply to the nucleon 
contribution, the only difference being that in the formula corresponding to Eq. ( BIO ), the factor e 2 L fJiU is replaced 
by 



CfL^v + KfMpv + 2efKfm f q 2 g f _ tL , . 



(B16) 



( f) 

By mimicking the steps that lead to Eq. (2.35), we obtain similarly looking formulas for lmn}f L , but with the functions 
Af, Bf, Df replaced by their imaginary parts, which in the end yield the results quoted in Section [VF. 
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